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Abstract 

Based on first principles calculations, we show that Mn impurities are magnetic in 
MgB 2 due to exchange-splitting of d 3z 2_ l band and they substantially modify B 
p a and p n bands through hybridization. Thus, Mn impurities could act as strong 
magnetic scattering centers leading to pair-breaking effects in MgB 2 . In contrast, 
we find Fe impurities in MgB2 to be nearly non-magnetic. 
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Magnetic impurities in a conventional phonon-mediated single band super- 
conductor [1] induce pair-breaking effects via spin- flip scattering [2]. However, 
for an unconventional multi-band superconductor, pair-breaking seems to be 
most probable from inter-band scattering [3,4,5], and thus it is not expected 
to depend solely on the exact atomic nature of the impurity. Since the dis- 
covery of superconductivity in MgB 2 [6], now confirmed to be a two-band 
superconductor [7,8,9,10,11,12,13,14,15], previous studies of chemical substi- 
tutions have been largely confined to examining the effects of inter-band and 
intra-band scattering of the electrons. However, due to the orthogonality of B 
p a and p n orbitals, the inter-band scattering effects which could lead to gap 
merging have not been observed in MgB 2 so far, although the superconducting 
transition temperature Tq showed a steady decrease as a function of increasing 
impurity concentration [5,16]. Of particular interests were Al [17,18,19,20,21] 
and 3d transition-metal substitutions [22,23,24,25,26,27,28,29] in the Mg sub- 
lattice and C substitutions [30,31,32,33,34] in the B sub-lattice of MgB 2 . 

Of all impurity substitutions studied so far, Mn shows the most rapid de- 
crease in T c as a function of its increasing concentration in MgB 2 [22,23,25]. 
However, no evidence of inter-band scattering was observed, instead the role 
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of magnetic exchange splitting of the bands in Mn-substituted MgB 2 was 
emphasized [23]. Earlier, it was shown that the inter-band scattering would 
be important if the lattice was distorted locally around the impurities [35]. In 
this regard, substituting C impurities for B atoms would have only a weak 
effect on inter-band scattering because it does not change the local symmetry. 
In contrast, substitutions in the Mg plane could create out-of-plane distor- 
tions and change the local point symmetry of nearby B atoms, resulting in 
a mixing of the in-plane p x ( y ) orbitals and the out-of-plane p z orbital. Under 
these circumstances significant a-n scattering would occur. 

In this letter, we show that Mn impurities are magnetic in MgB 2 due to 
exchange-splitting of d 3z 2_i band and they substantially modify B p a and p n 
bands through hybridization. In contrast, we find Fe impurities in MgB 2 to 
be nearly non-magnetic. Thus, the Mn impurities could act as strong mag- 
netic scattering centers leading to pair-breaking effects in MgB 2 . This may 
explain the rapid decrease of T c in Mg 1 _ x Mn x B 2 with increasing x, while 
in Mgi_ x Fe x B 2 a relatively slow decrease in T c with x may only be due to 
disorder-induced scattering of the electrons. Before we discuss our results in 
detail, we briefly describe the computational details. 

The normal metal electronic structure of the disordered alloys Mg x _ x Mn x B 2 
and Mgi^ x Fe x B 2 are calculated using the Korringa-Kohn-Rostoker (KKR) 
Green's function method formulated in the atomic sphere approximation (ASA) 
[36], which has been corrected by the use of both the muffin-tin correction for 
the Madelung energy, needed for obtaining an accurate description of ground 
state properties in the ASA [37], and the multipole moment correction to the 
Madelung potential and energy which significantly improves the accuracy by 
taking into consideration the non-spherical part of polarization effects [38,39]. 
Chemical disorder was taken into account by means of coherent-potential ap- 
proximation (CPA) [40]. The exchange and correlation were included within 
the local density approximation (LDA) using the Perdew-Wang parameteri- 
zation [41] of the many-body calculations of Ceperley and Alder [42]. During 
the self-consistent procedure the reciprocal space integrals were calculated by 
means of 1771 k-points in the irreducible part of the hexagonal Brillouin zone 
(BZ), while the energy integrals were evaluated on a semicircular contour in 
the complex energy plane using 20 energy points distributed in such a way 
that the sampling near the Fermi energy was increased. When calculating the 
density of states the number of k-points were increased to 2299. The Green's 
function was evaluated for 1000 points on a line in the complex energy plane 
parallel to the real axis and was then analytically continued towards the real 
axis. The atomic sphere radii of Mg and B were kept as 1.294 and 0.747, re- 
spectively of the Wigner-Seitz radius. The sphere radii of the 3d substituents 
were kept the same as that of Mg itself. The overlap volume resulting from the 
blow up of the muffin-tin spheres was approximately 10%, which is reasonable 
within the accuracy of the approximation [43]. The charge self-consistency it- 
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Figure 1. The spin polarized, total and sub-lattice resolved partial densities of states 
of Mgo.95Mrio.05B2 (top panel) and Mgo.95Feo.05B2 (bottom panel) alloys calculated 
as described in the text. The transition-metal d and B p densities of states for the 
respective alloys are also shown in the side panels. The dashed vertical line through 
energy zero represents the Fermi energy 

erations were accelerated using the Broyden's mixing scheme [44], the CPA 
self consistency was speeded up using the prescription of Abrikosov et al [45]. 
The calculations for Mgi- x Mn x B 2 and Mgi_ x Fe x B 2 alloys with x = 0.05 
were carried out at the experimental lattice parameters of MgB 2 [6], assum- 
ing a rigid underlying lattice. The effects of local lattice relaxation as well as 
any possible short-range ordering effects are not considered. Our results for 
Mgi- x Mn x B2 and Mg\- x Fe x B2 for x = 0.05 are analyzed in terms of to- 
tal, sub-lattice resolved partial densities of states, and Bloch spectral density 
A B (E, k) [46] evaluated at high symmetry points of the hexagonal BZ. The 
electronic structure of body-centered cubic (bcc) Fe is calculated with a sim- 
ilar setup using the experimental lattice constant a = 5.425 a.u.. The density 
of states of bcc Fe is calculated with 2470 k-points in the irreducible BZ. 

The calculated density of states (DOS) at the Fermi energy, Ep, of MgB 2 
was determined to be 3.493 states/Ry atom, which is consistent with previ- 
ous reports [47,48]. For Mgo.95Mno.05B2 and Mgo.95Feo.05B2 alloys, the spin 
polarized calculations showed a lower total energy when compared to the 
spin unpolarized calculations, suggesting the role of exchange splitting of the 
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Figure 2. The spin polarized, transition-metal ci-projected DOS in M g 9 ^MriQ^B 2 
(top panel) and Mgo.95-Feo.o5-B2 (bottom panel) alloys, resolved into three 
non-degenerate xy (x 2 — y 2 ), yz (zx), and 3z 2 — 1 symmetry components. The 
inset in the bottom panel shows the spin polarized, d-projected DOS of bcc Fe, 
resolved into two non-degenerate xy (yz,zx) (dashed line) and 3z 2 — 1 (x 2 — y 2 ) 
(solid line) symmetry components. In the inset figure, the x-axis corresponds to en- 
ergy ranging from —0.4 Ry to 0.2 Ry, while the y-axis represents the DOS from 
to 12 state/Ry/atom for both spins. The dashed vertical line through energy zero 
represents the Fermi energy. Note that the DOS in the figure corresponds to per 
transition-metal atom. 

states, in particular for Mn substituted MgB 2 . The magnetic energies for 
Mg 095 Mn 00 5B 2 and Mg 095 Fe om B 2 alloys are calculated to be 5.319 mRy 
and 0.001 mRy, respectively. The corresponding local magnetic moment for 
Mn is equal to 1.84 // B /atom, while Fe tends to be feebly magnetic with a 
local magnetic moment of 0.04 jU^/atom. In Fig. 1 we show the spin polarized, 
total and sub-lattice resolved partial densities of states of Mg 0i95 Mn 9 ^B 2 and 
Mg , 95 Fe om B 2 alloys. In the same figure we also show the transition-metal 
d and B p x ^ and p z DOS in these alloys. It is clear from Fig. 1 that the 
exchange-split Mn d orbitals in Mg 095 Mn , 05 B 2 alloys make Mn impurities 
magnetic in MgB 2 . We also find a substantial distortion in the majority and 
minority B p x ( y ) and p z bands due to their hybridization with the Mn d or- 
bitals. In addition, we find the B p band to have a very small exchange-splitting 
of 0.7 mRy in Mg 095 Mn om B 2 but no splitting in Mg 095 Fe om B 2 . 

To gain further insight into such a contrasting behavior of Mn and Fe im- 
purities in MgB 2 , we show in Fig. 2 the transition-metal ci-projected DOS in 
Mg $ 5 Mno.o5B 2 and Mgo. 95 Feo.o5B 2 alloys, resolved into three non-degenerate 
xy (x 2 — y 2 ), yz (zx), and 3z 2 — 1 symmetry components. We find that the 
band with 3z 2 — 1 symmetry is the narrowest of all the d bands. Thus, a 
half-filled 3z 2 — 1 band is more likely to induce exchange splitting than the 
other bands. In the case of Mn the minority 3z 2 — 1 band is nearly half-filled, 



4 



resulting in an exchange-split d band. Our calculations show that the centers 
of the d | and the d { bands are separated by 0.114 Ry leading to the Stoner 
exchange integral / (= ) to be 0.062 Ry//z.B. However, in the case of Fe 
the additional one electron fills up the minority 3z 2 — 1 band, as can be seen 
from Fig. 2, thereby making Fe nearly non-magnetic in MgB 2 . In addition, 
following the prescription of Liechtenstein et al [49], we have calculated the 
on-site Heisenberg exchange interaction J . The values of J for Mn and Fe 
alloys are found to be 0.7292 mRy and 0.0001 mRy respectively, confirming 
Mn to be magnetic and Fe to be nearly non-magnetic in MgB 2 . 

It is instructive to compare the DOS of Fe in Mg 095 Fe ,o 5 B 2 with that of 
bcc Fe as shown in the inset of Fig. 2, where the ci-projected DOS of bcc Fe 
has been resolved into two non-degenerate xy (yz, zx) and 3z 2 — 1 (x 2 — y 2 ) 
symmetry components. The calculated magnetic moment of bcc Fe is found 
to be 2.24 /is/atom. The most significant changes in the DOS of Fe involve 
x 2 —y 2 and 3z 2 — l symmetry components. Due to the change in local symmetry, 
as one goes from bcc Fe to hexagonal Mg 095 Fe , 05 B 2 , the narrow x 2 — y 2 
band of bcc Fe hybridizes with the B p x ( y ) bands of Mg , 95 FeomB 2 and gets 
broadened as shown in Fig. 2. The effects of hybridization on the B p x ( y ) DOS 
can be seen in Fig. 1 in the form of a hump around 0.14 Ry below the Fermi 
energy. On the other hand, the filling up of the minority 3z 2 — 1 band by the 
extra electron results in an inward movement of the band by 0.15 Ry, making 
magnetism unsustainable in Mg 095 FeQ^B 2 . 

Many of the superconducting properties of MgB 2 can be directly related to 
its normal state electronic structure. In particular, an understanding of how 
the additions of Mn and Fe impurities affect the electronic states over the 
entire BZ can be very useful in this regard. In Fig. 3 we show the spin po- 
larized, k-resolved Bloch spectral density Ab(E, k) for Mg Q95 Mn x) 5 B 2 and 
Mg 0i95 FeomB 2 alloys at T, A, M, H, K, and L points of the hexagonal BZ. As 
expected, the differences between majority and minority Ab(E, lc) are much 
more for Mg 095 Mn _ 05 B 2 alloys than for Mg 095 Fe 005 B 2 alloys over the entire 
BZ. New states, predominantly d in character are seen to form just below E F 
at T and K points. The disorder-induced broadening of the doubly degenerate 
B p a states located between —0.05 < E < +0.05 at T and B p n states just 
above Ep at M point can be seen clearly. At M point, the lowest lying peak 
is characteristic of the s band, while the next two peaks correspond to the p a 
bands. We find that the p a bands at M point in ordered MgB 2 are relatively 
far apart in comparison to the Mn and Fe substituted systems. A change in 
the dispersion of bands is expected to affect the three dimensionality of the 
Fermi surface [47] of the substituted materials (a detailed comparison will be 
published elsewhere). We also note that in transition-metal diborides, such 
as TaB 2 and others, a strong metal d hybridization with the B p bands is 
observed in A-L-H plane of the BZ [50,51]. However, in the case of Mn and 
Fe impurities, the d bands are narrower and only the states that are close to 
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Figure 3. The spin polarized, k-resolved Bloch spectral density Ab(E,\s) for 
M go.95Mno.05B2 (left) and Mgo.95Feo.05B2 (right) alloys at T, A and M (top panel), 
and H , K, and L points of the hexagonal Brillouin zone. The dashed vertical line 
through energy zero represents the Fermi energy 

Ep are affected. 

In summary, we have carried out density-functional-based electronic structure 
calculations to study the effects of Mn and Fe substitutions on the a and 
7r bands of MgB 2 . The self-consistent calculations for Mg _g 5 Mn _ 05 B 2 and 
Mg , 95 Feo,Q 5 B 2 alloys show that Mn forms a local magnetic moment of 1.84 
/is/atom while Fe tends to remain feebly magnetic with a local magnetic 
moment of 0.04 ^s/atom. Further, the k-resolved Bloch spectral density of 
states at T, A, M, H, K, and L points of the hexagonal Brillouin zone show 
that for Mn substituted alloys the a bands of MgB 2 are strongly hybridized 
with the Mn d bands and undergo considerable spin-dependent modifications 
when compared to its Fe counterparts. Given that p a and p n orbitals are 
responsible for superconductivity in MgB 2 , Mn impurities could, therefore, 
act as strong magnetic scattering centers leading to pair-breaking effects. In 
contrast, the Fe impurities in MgB 2 are nearly non-magnetic with essentially 
no pair-breaking effects. 

One of us (PJTJ) would like to thank Dr. Andrei Ruban for providing the 
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KKR-ASA code. Discussions with Dr. Igor Mazin on the electronic structure 
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